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Cosmic Microwave Background Spectral Distortions from Cosmic String Loops 
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Cosmic string loops contain cusps which decay by emitting bursts of particles. A significant 
fraction of the released energy is in the form of photons. These photons are injected non-thermally 
and can hence cause spectral distortions of the Cosmic Microwave Background (CMB). Under the 
assumption that cusps are robust against gravitational back-reaction, we compute the fractional 
energy density released as photons in the redshift interval where such non-thermal photon injection 
causes CMB spectral distortions. Whereas current constraints on such spectral distortions are not 
strong enough to constrain the string tension, future missions such as the PIXIE experiment will be 
able to provide limits which rule out a range of string tensions between Gfi ~ and G/r ~ 

thus ruling out particle physics models yielding these kind of intermediate-scale cosmic strings. 


I. INTRODUCTION 

Cosmic strings [1] are topologically stable solutions of 
the classical field equations in many particle physics mod¬ 
els beyond the Standard Model. If matter is described by 
such a model, then a simple causality argument of Kib¬ 
ble [2] ensures that a network of strings will form during 
a symmetry-breaking phase transition in the early uni¬ 
verse and will persist to the present time. Both in the 
1980s (see e.g. [3-5] for reviews) and more recently (see 
e.g. [6] for a recent review) there has been a lot of work 
devoted to the observational signatures of such strings. 
The network of cosmic strings is described by a single 
free parameter, namely the mass per unit length /r of 
a string, which is related to the energy scale rj of the 
symmetry-breaking phase transition via ^ = drj^, where 
d is a number of order unity (which, in particular, to 
a first approximation does not depend on coupling con¬ 
stants in the particle physics model). 

Cosmic strings carry energy and hence their gravita¬ 
tional effects can lead to cosmological signatures. In fact, 
at one point, it was conjectured [7] that cosmic strings 
could seed all of the observed structure in the Universe. 
With the discovery of acoustic oscillations in the angular 
power spectrum of the Cosmic Microwave Background 
(CMB) [8] this possibility was ruled out. Currently, the 
best and most solid constraints on the string tension 
come from precision measurements of the CMB angular 
power spectrum. They are [9, 10] 

Gfi < 2x 10”'^ (1) 

(see [11] for older results), where G is Newton’s gravita¬ 
tional constant (the string tension is typically expressed 
in terms of a dimensionless number by multiplying by 
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G).^ This already rules out a set of particle physics mod¬ 
els which have a symmetry breaking scale at the high end 
of the “Grand Unification” class. Tightening these con¬ 
straints or providing new constraints is of great interest 
since it will provide new ways to test particle physics 
beyond the Standard Model, ways which are comple¬ 
mentary to the test at accelerators such as the Large 
Hadron Collider (LHC) (see [12] for an elaboration on 
this theme). 

Cosmic strings also form closed string loops. Kibble’s 
causality arguments coupled to studies of the scaling so¬ 
lutions of the Boltzmann-type equation which describes 
the time evolution of the network of strings show [3-5] 
that there will at all times be a scaling network of infinite 
strings with typical separation and correlation length of 
the order t (the cosmic time), as well as a distribution 
of string loops which are the result of interactions of the 
long strings. This distribution of string loops is statisti¬ 
cally independent of time if all lengths are scaled to the 
Hubble radius t. 

There has been a lot of work on the gravitational ef¬ 
fects of cosmic strings (see e.g. [6] for a recent review). 
The long strings lead to distinctive line discontinuities 
in CMB temperature maps [15], to patches in the sky 
with direct B-mode polarization [16] and to wedges in 
21cm redshift maps at high redshifts with extra absorp¬ 
tion [17]. Cosmic string loops may form the seeds of 
globular clusters [18], they may assist the formation of 
high redshift super-massive black holes [19], and they 
may lead to ultra-compact mini-halos embedded within 
galaxies [20]. All of these effects are due to the gravity 
of cosmic strings. 

In this note we study a non-gravitational effect of cos¬ 
mic string loops. As argued in [21] and recently con¬ 
firmed numerically in [22], cosmic string loops are not 
smooth but generically contain cusps. Cusps are unsta¬ 
ble to decay and will evaporate into jets of particles. A 
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significant fraction of the energy lost in this process goes 
into non-thermal photons. Non-thermal photons gener¬ 
ated during in the redshift range 3 x 10® > z > 3.6 x 10® 
(where z denotes cosmological redshift) are not able to 
reach equilibrium with the CMB. As a result, these pho¬ 
tons cause spectral distortions of the CMB which can be 
measured by future experiments. We show below that, at 
least for a class of cosmic string models, results from the 
planned PIXIE experiment could rule out string loops 
with string tensions in the range 10“®® < Gfji < 10®^. 

Consequences of cusp annihilation for cosmology have 
been studied earlier in the context of ultra-high energy 
cosmic ray and neutrino production from strings [23]. 
These works did not consider spectral distortions. Spec¬ 
tral distortions caused by strings were considered in [24] 
in the context of decaying cosmic defects, in [25] in the 
context of acoustic waves produced by cosmic strings, 
and in [26] in the context of superconducting cosmic 
strings. 

Effects of the final collapse of a string loop (once the 
radius becomes comparable to the width) on black hole 
formation [27] and on dark matter production [28] have 
also been studied. Note that the photon production dur¬ 
ing the final loop collapse is negligible in magnitude com¬ 
pared to the amount of energy released from string cusps, 
and hence we do not consider this final collapse in this 
paper. 


II. ENERGY LOSS FROM CUSP 
ANNIHILATION AND GRAVITATIONAL 
RADIATION 


We will work in the context of a simple one-scale model 
for the cosmic string loop distribution according to which 
loops are formed at times t with a radius 

= f, ( 2 ) 

where a and (3 are constants. The average length of a 
string loop is I = PR, and the constant a relates the 
string length at the time of formation to the time. Ana¬ 
lytic arguments suggest /3 ~ 10, while numerical cosmic 
string evolution simulations [13] give the rough estimate 
a ~ 0.1. Within a radiation-dominated epoch the re¬ 
sulting number density of string loops per unit radius is 
given by 


n{R,t) = Va®/ 2 ; 3 “®/ 2 i- 3 / 2 ^- 5 / 2 ^ 

for Rc{t) < R < (3) 

Here the constant N is determined by the number of long 
string segments for Hubble volume. Numerical simula¬ 
tions indicate that this constant is of the order N ^ 10 
[13]. The origin of the lower bound Rc{t) will be ex¬ 
plained momentarily. 


Once formed, loops slowly decay through numerous 
channels,^ and the number density (3) does not include 
these effects. As a simple way to incorporate these decay 
processes into n{R,t), we impose a lower cutoff Rc{t) on 
the radius in the distribution of string loops. This cutoff 
is the radius for which a string loop will evaporate com¬ 
pletely within one Hubble time scale. Below the cutoff 
the number density is roughly: 


Within our model there are two competing decay mech¬ 
anisms: decay through the emission of gravitational ra¬ 
diation [14], and decay through cusp annihilation. The 
dominant mechanism determines the cutoff scale Rc{t). 

The emission of gravitational radiation is independent 
of time as well as the radius of the loop and satisfies the 
simple relation 

R = -yG^, (5) 


where the dot denotes a derivative with respect to t 
and where 7 is another constant which is related to the 
strength of gravitational radiation and which must be de¬ 
termined from numerical simulations. Its value is of the 
order 7 ^ 10. Noting that the energy contained in a loop 
is E = = fJ,pR, the power emitted by a loop due to 

gravitational waves is thus 


dE 

dt 


GW 


PlG^l^. 


( 6 ) 


The critical radius may then be estimated by Rc = —Rt 
which yields 

R!p^{t) =-fG^xt. (7) 


Next we analyze loop decay through cusp annihilation. 
When speaking about string cusps, one is working in the 
Nambu-Goto approximation of cosmic string dynamics in 
which the hnite thickness of the string is neglected. This 
is an excellent approximation in the context of cosmol¬ 
ogy since the spatial extent of the string is cosmological 
whereas the width is microphysical. However, at a string 
cusp the Nambu-Goto approximation breaks down. If we 
consider a string loop with characteristic length scale R, 
then, as shown in [29], the two string segments of the 
loop around the cusp will overlap (i.e. their distance will 
be smaller than the string width) over a distance 

1{R) = kR^/^wG^ , (8) 

where n is an order unity constant and w is the string 
width. This width is related to the string tension by 
w ^ ? 7 “® ^ /i“®/^. This region contains an energy of 

6E = 2nl{R) = 2KnR^GyjG3 _ ( 9 ) 


^ The defects we consider in this paper are topologically stable. 
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Locally, a cusp looks like an overlapping string and an- 
tistring segment. No topology protects this region of the 
string from exploding into a burst of particles. The pri¬ 
mary particles in this burst will consist of quanta of the 
scalar and gauge fields which make up the string. The 
corresponding particles are unstable and will decay into 
a jet of stable Standard Model particles. We call this the 
cusp annihilation process [30]. Most of the energy will 
end up in pions, neutrinos and photons. We expect that 
a fraction / of order one will end up in photons. Since 
the strings are out-of-equilibrium field configurations, the 
photons resulting from cusp annihilation will be very en¬ 
ergetic and out of thermal equilibrium. Photons emitted 
in the redshift range 3 x 10® > z > 3.6 x 10® will not be 
able to thermalize with the CMB and will hence lead to 
CMB spectral distortions. Observational limits of CMB 
spectral distortions hence can lead to limits on the cosmic 
string model. 

Since there is of the order one cusp per loop oscillation 
time, and since the loop oscillation time is of the order 
of the radius R, the rate of energy loss of a string loop 
due to cusp annihilation is 


AE 

"Af 


cusp 


2 K/i ^ 


w\ 1/3 

r) 


( 10 ) 


It follows that the rate of change in the loop radius due 
to cusp annihilation is 


R = - 


2k 

T 



( 11 ) 


The decay timescale for a loop of initial radius Ri may 
be obtained from the approximation 

R{t) ~ tR{ti) -\- Ri. (12) 

The resulting cutoff in the loop distribution is 

(13) 

We may now compare the cutoff radii due to gravi¬ 
tational radiation and cusp annihilation. These cutoffs 
scale with time like R^^ {t) ~ t and i?™"P(t) ~ 
respectively. Thus, at early times cusp annihilation pro¬ 
vides the larger cutoff, while at sufficiently late times the 
cutoff due to gravitational radiation is larger. The cutoffs 
are equal at the “cross-over” time T 


8k®w _ 8k®G1/2 1 

/ 3 ® 7 V^G 4 “ /3374AI/2 (^G) 9/2 ■ 


(14) 


In the last equality we have made the relationship be¬ 
tween w and p precise by inserting w = 
where A is a model-dependent Higgs coupling constant 
which depends on, e.g., the coupling constant of the QFT 
generating the cosmic strings. The cross-over time may 
also be described in terms of the power emission as fol¬ 
lows: the power emitted through gravitational radiation 


and cusp annihilation, i.e. (6) and (10), become equal at 
a radius 


8k®G1/2 

7®/33Ai/2(/rG)7/2- 


(15) 


For R > Ri, gravitational radiation is more efficient, 
while for R < Ri, cusp annihilation is more efficient. The 
time T is the time when R'^'^{T) = i?™®P(T) = i?*. 

To summarize, we take as our cutoff of the loop number 
density distribution n{R, t) 

Rc{t) = i?™"P(t)0(r -t)+ R^'^[t)Q{t - T) 

= 7 mG -t) + - T)^ , (16) 


where as usual 0(a;) denotes the Heaviside function. As 
we will see in the next section, the energy density in 
photons produced from cusp annihilation is dominated 
by loops with radius near the cutoff radius. 


III. CALCULATION OF PHOTON 
PRODUCTION FROM CUSPS 


We now turn to the calculation of the photon en¬ 
ergy density input from cusp annihilation. We assume 
that a fraction / of the cusp energy ends up in non- 
thermal photons. Photons released at redshifts larger 
than Zin ~ 3 X 10® are able to thermalize. Hence, we are 
interested in redshifts between Zin and the the redshift of 
re-ionization. To slightly simplify the algebra, we will in 
fact focus on redshifts larger than the redshift of matter 
radiation equality which is Zgg = 3.6 x 10®. The corre¬ 
sponding times are ti„ and t^q. We will focus on some 
time t in the window tin < t < teq- 

All cusps present between tin and t will contribute to 
the non-thermal energy in photons. At each time t' be¬ 
tween tin and t, all loops present at that time will radi¬ 
ate photons via cusp annihilations. The energy density 
in photons produced at time t' redshifts like radiation 
between t' and t. Since there is roughly one cusp for 
loop oscillation time R, the non-thermal photon energy 
density is 

rat'/P i 

X / dRn{R,t')—f2KpR^^^w^^^. 

Jo R 

Recall that our calculations are in the radiation- 
dominated phase, so a{t) oc It will be convenient 

later to have this quantity normalized with respect to 
the background (photon) energy density Pb{t). From the 
Friedmann equation we easily compute 


Pb{t) 


_3_ l_ 

■ 


(18) 
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Combining these equations we have 


Pb{t) 


647r 




Jo 


at' IP 


dRn{R, t')R 


- 1/3 


(19) 


The integral over R may be approximated by 

at' IP 

dRn{R,t')R-^^^ 

/ \ 5/2 „ 

»N ( 1 ) ( 20 ) 

We obtain this by considering the integral over the range 
Rc{t) < R < at/fJ. In this region the loop distribution 
n{R,t) behaves like R~^l'^, and as a result the integral 
is peaked near R ~ Rc(t). There is also a contribution 
from the integral over the region 0 < R < Rc{t) which 
is, at most, of the same order of magnitude. Since the 
precise form of the loop distribution in this region is not 
known, we drop this contribution with the understanding 
that our calculation is accurate only to within an order 
of magnitude. Inserting (20) into (19) and tidying up we 
obtain 


/ dt'{t')^/^Rc{t')-^^/^ , ( 21 ) 
Pb[t) Jtir, 


where we have defined the constant 


Cl 


12871 

11 



5l2 

TV/kA-Cs. 


( 22 ) 


The remaining integral over t' may be performed ex¬ 
actly. For our purposes we are interested in maximiz¬ 
ing the fraction pj(t)/pb{t). Since the integrand in (21) 
is positive, we achieve the maximum value by letting 
t —>■ teq. It is also convenient to recast factors of the 
string tension p,G in terms of the cross-over time T (re¬ 
call T oc [plG)~^^^). Inserting our expression (16) for the 
cutoff radius and tidying up further we obtain 


Pyjteq) 

Pbiieq) 


G^G^IWr'^Io 



dt' 




t') 


+ -T)], (23) 


where the overall constant is now 




- 2/9 


(24) 


Case (i): T < tin- For this case Rc{t) = R^^{t) 
over the entire range of integration. After integrating we 
obtain 


P'ljteq) 

Pb{teq) 


3G2GC18t2/9 


{tin) - {teq) 


(25) 


Note that this scales like and so is maximized by 

letting T tin- 

Case (ii): t^q < T. For this case Rc{t) = i?™®P(t) 
over the entire range of integration. Integrating yields 
the expression 


P'fjteq) 

Pb{teq) 


8G2GCI8T-IC72 


(Gg)C 8 - (ti„)C 8 


(26) 


This expression scales like T and so is maximized 

by letting T -)■ 4,. 

Case (Hi): Cn < T < tf,q. For this case the functional 
form of the cutoff radius Rc{t) changes within the inte¬ 
gration region. The result of integration is 


P'fjteq) _ ^ 
Pbiteq) ~ " 


rp-1/0 



+ 8T-1C24 


T^/^-{Un)^/^]]. 

(27) 


This regime of T contains the configuration which max¬ 
imizes p^{tf.q)/Pb{teq)- We denote by Tmax the value of 
T which maximizes p^{teq)/Pb{teq)] this value is found 
numerically to be 

^ = 6.6 X 10^°. (28) 

VG 

Note that this value is determined by tin and teq alone; it 
does not depend on any other parameters of our model, 
all of which have been sequestered into the constant C 2 ■ 
In terms of redshift, T^ax translates to Zmax = 7.4 x 10®; 
from this we see that T^ax is quite close to the initial 
time. The maximum value of the non-thermal photon 
energy density is thus 

= G2(1.3 X 10"®). (29) 

Tinax 

We conclude this section by discussing the numerical 
value of G 2 . This constant is a conglomeration of the 
many constants of our model; explicitly, G 2 may be writ¬ 
ten 


Pjjteq) 

Pb{teq) 


G 2 ~ 23 iV/a®/2/3-ii/6^-i7/i8^i/3;^-i/i8_ ( 30 ) 

In the text above we have given the order-of-magnitude 
estimates for these parameters as 

iV = /3 = 7=10, a = / = 0.1, A = 7v = 1, (31) 


There are now three possible cases, depending on whether which yields the value 
the cross-over time T lies below, within, or above the 

range of integration. G 2 = 1.2 x 10“"^^^ 


( 32 ) 
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We include in C 2 an uncertainty of two orders of mag¬ 
nitude. This is due mainly to the uncertainty in the 
parameters a and JV which need to be determined by nu¬ 
merical cosmic string simulations. Combining (29) and 
(32) we obtain the approximate maximum value of the 
photon density 

= 1.6 X 10-9±2. (33) 

Pb[teq) max 

The value of the string tension associated to this maximal 
value is 

^GUax ~ 2.2 X 10-^^ (34) 

In the next section we turn to the implications of these 
results for the magnitude of CMB spectral distortions. 

IV. CURRENT AND FUTURE CONSTRAINTS 
ON THE COSMIC STRING TENSION 


The proposed PIXIE experiment [34] will improve the 
constraint by four orders of magnitude to 

^ < 10-9. (40) 

Pb{t) 

In the previous section we found that for the optimal 
value of T the injected photon density is 1.6 x lO”®^^. 
This value is too small to be constrained by COBE data; 
however, it does potentially lie within the sensitivity of 
PIXIE. Taking a large value for C 2 (i.e., considering the 
+2 in (32)), our model predicts that PIXIE would detect 
CMB distortions due to cosmic string loops with string 
tension in the range 

2.2 X < fiG < 9.4 X 10"^^ (41) 

Smaller values of C 2 result in a smaller range, and suf¬ 
ficiently small values of C 2 (e.g., considering the —2 in 
(32)) result in CMB distortion levels below the sensitivity 
of PIXIE. 


Non-thermal injection of photons in the redshift in¬ 
terval between z ~ 3 x 10® and recombination leads to 
spectral distortions in the CMB (see e.g. [31] for a recent 
comprehensive overview, and [32] for an earlier paper). 
For 3 X 10® > z > 10®, the energy injection leads to a 
Bose-Einstein distribution modified by a chemical poten¬ 
tial (^ distortion). The chemical potential /i generated 
by energy injection is given by 


1 SU 

~ ^ir 


(35) 


where SU is the energy density in injected photons, and U 
is the background photon density. In our case, SU = p-y, 
the energy density from cusp evaporation computed in 
the previous section, and U = pb- 

For redshifts in the range 10® > z > 3.6x 10® (the lower 
bound being the redshift of matter-radiation equality), 
the energy injection produces a Comptonized spectrum 
characterized by a y-distortion, with the y-parameter 
given by 


iir ■ 


(36) 


The best limits on spectral distortions still come from 
the COBE experiment, and are [33] 

1^1 <9x10-®, (37) 


and 


|y|< 15x10-®, (38) 


respectively, both at 95% confidence level. These bounds 
give essentially the same order of magnitude constraint 
on the injected photon energy density, which is thus in¬ 
dependent of time: 

< 10-®. (39) 


V. CONCLUSIONS AND DISCUSSION 

Cusp annihilation is a mechanism by which cosmic 
strings lose energy to photons. Within our model of 
cosmic string evolution, cusp annihilation is the domi¬ 
nant decay channel at early times. After the cross-over 
time T, cusp annihilation gives way to gravitational radi¬ 
ation as the dominant decay channel. The photons pro¬ 
duced by cusp annihilation are out-of-equilibrium and 
can lead to spectral distortions of the CMB. We have 
computed the energy density of this injected photon flux 
as a function of the cross-over time T (equivalently, as 
a function of the string tension pG). We find that the 
resulting CMB distortion is largest for a cross-over time 
T at redshift ~ 7 x 10®, corresponding to a string tension 
/rG - 2 X 10-1®. 

Current limits on p and y distortions of the CMB are 
not sensitive enough to constrain cosmic string models. 
However, our analysis shows that the planned PIXIE mis¬ 
sion will have the sensitivity to rule out (or confirm) 
string tensions in the range of 10-®® < Gp < lO-^^ 
which corresponds to symmetry breaking scales 77 of 
10“GeV <r]< 10i®GeV. 

An important caveat is that our analysis does not take 
the gravitational back-reaction of cusps into account. 
Back-reaction may prevent cusps of the length predicted 
by the Nambu-Goto approximation to develop. In this 
case, less photons would be produced and the constraints 
on the cosmic string tension would be weaker. 

Note added: After the initial submission of this work 
we became aware that spectral distortions from cosmic 
string loop cusp decay were already considered in Section 
4.2 of [35] (which is based in part on [36]). Our qualita¬ 
tive results agree, although the quantitative details differ 
a bit due to different assumptions made about the struc¬ 
ture of a cusp. This difference highlights the importance 
of improved studies of back-reaction effects on cusps 
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